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Abstract 

Following the Caldeira-Leggett approach to describe dissipative quantum systems the structure 
function for a harmonic oscillator with Ohmic dissipation is evaluated by an analytic continuation 
from euclidean to real time. The analytic properties of the Fourier transform of the structure 
function with respect to the energy transfer (the "characteristic function") are studied and utilized. 
In the one-parameter model of Ohmic dissipation we show explicitly that the broadening of excited 
states increases with the state number without violating sum rules. Analytic and numerical results 
suggest that this is a phenomenologically relevant, consistent model to include the coupling of a 
single (sub-)nuclear particle to unobserved and complex degrees of freedom. 
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I. INTRODUCTION 



The structure functions measured in deep inelastic lepton scattering provide important 
information about the nature and momentum distribution of the constituents of the target. 
This is particularly true in the relativistic domain where the point-like building blocks of 
matter, the quarks and gluons, have been discovered and studied by inclusive scattering of 
multi-GeV electrons, muons or neutrinos. In a non-relativistic description the (longitudinal) 
structure function for one scalar particle is given by 

S(q,u) = £ 6(v-(E n -E )) |<n|exp(zg-x)|0>| 2 (1) 

where the summation is over all discrete and continuum states which are excited by the 
probe. Here q, v denote momentum and energy tranfer and E n the excitation energies of the 
target. Very often the confinement of quarks is described by using rising potentials which 
lead to a purely discrete spectrum. The simplest version is given by a harmonic oscillator 
potential which has a structure function 

The oscillator frequency and length are denoted by ojq and b = 1 / y/mtuo where m is the 
mass of the struck particle. For simplicity, all many-body and recoil effects have been 
neglected and the elastic line (n = 0) is included in the sum over excited states. However, the 
observed structure functions are smooth due to hadronization and/or final state interactions. 
A number of recent theoretical studies have accounted for that by simply smearing out the 
(^-functions in eq.(J2J) by a Breit-Wigner distribution with a constant width or averaged over 
nearby u-bms |2|. It is obvious, that this is not only ad hoc but also may violate general 
properties of the structure function, e.g. the fact that it has to vanish below the first 
(positive) excitation energy. 

It is the purpose of this paper to demonstrate that a consistent quantum mechanical 
framework exists which allows to treat couplings to unobserved degrees of freedom in a 
simple manner. Several methods have been used in the past to achieve that, for example for 
Coulomb excitation of particle-unstable states jj. For quasielastic scattering of electrons 
from nuclei Horikawa et al. ^ first have included multi-nucleon channels by employing an 
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optical potential without violating the non-energy-weighted sum rule (NEWSR) 

POO 

/ dvS(q,v) = 1. (3) 

J 

However, there exists a simpler treatment based on the description of dissipative quantum 
systems within the path integral formalism jsj]. This originates in the celebrated work of 
Feynman and Vernon (| and Caldeira and Leggett [tJ who have modelled the coupling of 
the system to an environment of N(—>- oo) harmonic oscillators 



2 N / 2 -i \ N N 2 

H = £ + ^) + E U^ + 2 m ^ )- X -^ CnXn + x2 ^2^ (4) 

*ll<> n=1 \ tlUn * J n= i n= i ^" l >n UJ n 

with a bilinear coupling. The limit N — > oo of the number of environmental oscillators is 
essential in preventing the bounded motion of all particles to come back to the initial state 
after some time : the so-called Poincare recurrence time then tends to infinity [8] and irre- 
versibility becomes possible - still in a unitary quantum-mechanical framework. The infinite 
number of degrees of freedom also allows for strong damping even if each environmental 
oscillator couples only weakly to the system. This mechanism leads to a broadening (and a 
shift) of the (^-functions in the structure function of the confined system without violating 
the sum rules. 

In this approach the path integral description of the system offers particular advantages 
since the bath oscillators can be integrated out exactly. This gives rise to a retarded two- 
time action for the single particle which does not have a Hamiltonian counterpart anymore 
similar as in the time-honoured polaron problem j^]. In particular, no Schrodinger equation 
for the single particle motion is available. However, if this particle moves in a harmonic 
potential 

V(x) = -muo 2 x 2 (5) 

then the remaining path integral can also be done exactly. In section 2 we will employ this 
formalism and the explicit results for the damped harmonic oscillator to obtain the Fourier 
transform of the structure function with respect to the energy transfer. As thermal physics 
lives in euclidean times one needs an analytic continuation to real times which is performed 
in section 3. Numerical results are presented in section 4 while the conclusions are given in 
the final section. Some technical details are collected in two appendices. 
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II. CHARACTERISTIC FUNCTION AND CORRELATION FUNCTIONS OF 
DISSIPATIVE SYSTEMS 



We will calculate the structure function from its Fourier transform 



1 r+oo 

S(q,v) =: — / dte wt %{t), (6) 

111 J~oo 

the "characteristic function" Q]. For the pure harmonic oscillator one has 



\?hi (i 



— iuot 

Y ..... 



(7) 



which after expansion and integration leads to the result given in eq. (j2J)- For the damped 
harmonic oscillator the characteristic function can be related to the particular correlation 
function in euclidean time r 

T(q,r) = (o|T(e-^ (T) e^ (0) )|o) (8) 

by an analytic continuation. From the spectral representations 

*' W l = EI<»|exp(w i)|0>r( (9) 
T(q,T)\ n \ex V \-(E n -E a )\ T \} 

one sees that both expressions coincide for positive euclidean time and the correct analytic 
continuation is therefore obtained by considering T(q, r > 0) and replacing r — > it. 
In the path integral approach we have 

JVx exp[-zg ■ (x(r) - x(0))} exp{-^l[a:]} 
Kq ' } /J™ JVxexp{-A[x}} ' K } 

Here A[x] is the effective action of the particle after the oscillators of the environment have 
been integrated out and the limit (3 — > oo of the final euclidean time projects out the ground 
state of the system Q|. In this limit the boundary conditions for the path integrals in 
eq.(JTUJ) do not matter; therefore we may set x(—j3/2) = x(/3/2) = x and integrate over 
x, i.e. perform the thermodynamical trace. This allows us to directly take over results 
from dissipative quantum systems where similar correlation functions (e.g. for the position 
operator) have been evaluated at finite temperature, i.e. finite (3 Since this is quite 
standard now we can be brief and immediately use results from the nice review by Ingold 
jl^j . in particular from chapter 4.3 with the driving force F(a) = iq [5(a — r) — <5(er)]. In the 
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limit j3 — > oo the sum over Matsubara frequencies turns into an integral so that the final 
result for the euclidean correlation function reads 

q 2 i ^ +co i_ cos (£ r ) 



T(q,r) = exp 
Here 



dE 

2mnJ-oo E 2 + \E\-i(\E\)+ul 
J{uj) E 



(11) 



l{E) = / du^^— 12 

IK J nm Jo uj E + to v ' 

is the damping kernel which is produced by the coupling of the system to the heat bath. 
One does not have to specify all masses, frequencies and coupling constants in Eq.(jlJ) but 
only the spectral density J {(J) of the environment oscillators. The simplest assumption is 
Ohmic dissipation 

Johm = mjuj =$> 7ohm(£) = 7- (13) 

Although some observables (e.g. the ground state energy) need high-frequency cut-offs 
this form can be used with impunity for the structure function where only energy differences 
matter. For simplicity it will also be employed in the following. We then obtain 

T(q,r) = exp{-2q 2 [x v (0)-x v (r)}} (14) 

with 

x v (r) = — HdE— ^^1— , (rreal) . (15) 

VK ' 2m7rio E 2 + -fE + ui ' v ; v ; 

The analytic continuation of the above function which coincides with xv(t) for positive 
euclidean time 

Cv(r) = x v (t) , (r>0) (16) 



may be called the Vineyard function since many years ago Vineyard [12| derived an identical 
form of the structure function for inclusive scattering of slow neutrons from quantum liquids. 
Its nice feature is the clean separation between the squared momentum transfer and the 
variable t which is conjugate to the energy transfer. Such a form was also obtained in a 
(zeroth order) variational calculation of relativistic deep inelastic scattering from a scalar 

Hticle where the broadening of the elastic line was due to multiple meson production 
. We thus have found a simple description of inclusive scattering with one additional 
parameter 7 which accounts for the coupling of the particle to additional degrees of freedom 
|l8j . These could be the continuum and/or many-nucleon emission in the case of quasielastic 
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scattering from nuclei or the production of colorless hadrons in the case of scattering from 
quarks. 

Note that T(q, r = 0) = 1 so that the sum rule (jSJ) is fulfilled. We will also see that the 
resulting structure function has the correct support, i.e. no unphysical excitations occur. 
This is because the Caldeira-Leggett model is based on a consistent many-body Hamiltonian 
and the environmental degrees of freedom have been integrated out without approximation. 



2l\ or time 



In contrast, other descriptions of the damped harmonic oscillator [2Cj, friction 
asymmetry |22j in general require a modification of usual quantum mechanics. 

III. ANALYTIC CONTINUATION 



The main task left over is to perform the analytic continuation to get the "characteristic 
function" $ g (t) from the correlation function T(q, r). This requires the analytic continuation 
of the Vineyard function £y(r — > it). 

It should be emphasized that xy{r) and £y(r) are different functions which only coincide 
for r > 0. In particular, although from eq. (fTK|) xy(r) is even in r we will see that the 
Vineyard function has a logarithmic cut on the negative real r-axis. In the following we 
will always indicate the range of validity of the representation for the Vineyard function in 
parenthesis as was done in eq. (fT^|) . It is obvious that this form cannot be used for analytic 
continuation to Minkowski time (where scattering occurs) because the cosine-function would 
blow up. As shown in Appendix A we may, however, distort the integration contour as in 
and obtain for r > 



I I J in a .-Et / jt? „/ t?\ -Et 



ref. |2i 

^ W = 2^i« ^ iiq^F ^ =: L ™^ (n«r >o). (17) 

This now allows an analytic continuation r — > it to obtain the characteristic function 

$ g d - h -°-(t) = exp { -2q 2 [£ y (0) - fr(if) ] } • (18) 

To see the effects of the damping we insert eq. (fTTj) into eqs. (fT8|).([fi|) , expand the exponential 
and perform the t-integration. This gives 



S d - h - -(q,u) = e- qb/2 5(v) + Y — — e 

roo I n \ 

■J Q dE 1 . . . dE n 5\v~Y, E kj P( E i) ■ ■ ■ P( E n) ■ (19) 



The first term is the elastic line with the square of the typical gaussian form factor for the 
harmonic oscillator. However, the oscillator length is renormalized by the interaction with 
the environment: 

b 2 — & 2 = 4^(0) = b 2 ^-arctanf— ) < b 2 (20) 

\l n V 7 / 

with 



fi = ^L 2 -^. (21) 

(we are only considering the underdamped case 7 < 2uo ) • Compared with eq. (J2J) all excited 
states are now broadened; in particular, the (n = l)-term just maps the weight function 

2g 2 e-« 2fc2 / 2 Q(u) p{v) = ^ e^/ 2 -**>- ^ 7 ffiL - . (22) 

H V JHK ' 2 u + u 2tt (i/-cu ) 2 + r 2 (i/)/4 V ' 

Apart from the additional factor 2uq/ (v+uq) the line shape is just an one-sided Breit-Wigner 
distribution with the energy-dependent width 

r < £ > = i^- < 23 > 

Note that this distribution vanishes at threshold and is identically zero for unphysical neg- 
ative energy transfers. The same holds for all other terms in the expansion (jl9j) due to 
the 5-function and the fact that all Ek > 0. Thus the structure function of the damped 
harmonic oscillator has the correct support. Unfortunately, it is not possible to evaluate the 
higher-order terms analytically. Only in the narrow-width approximation 

T(E) « I>o) = 7 (24) 

a simple result is found when additionally the prefactor is also evaluted at v — ouq and the 
E- integration extended to —00: 

i Y {it) « -jexp [~iu t - ||tf) . (25) 

Note that an exact expression of this form together with a real correction r(t) is derived in 
Appendix |X] One then sees from eqs. (|THjl.(|5jl that the 5-function in the n-th excited state 
contribution to the structure function (J2J) of the harmonic oscillator turns into a Breit-Wigner 
function with a width 

r n « n 7 . (26) 
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This result is well-known from the density of states of the damped harmonic oscillator 24 1 



and has also been discussed for the width of multi-phonon giant resonances [25J. In the 
present context, however, it should be stressed that it is only approximate and leads to a 
small, but non- vanishing structure function for v < 26]. This is due to the wrong analytic 
behaviour of the Vineyard function in the approximation (J25)) where \t\ = also produces 
a cut for Im t < 0. In contrast, the exact expression has only a logarithmic cut in the upper 
half t-plane in accord with general properties of the characteristic function. This can be 
best seen in the explicit expression of the Vineyard function with Ohmic damping in terms 
of the standard exponential integral which is derived in Appendix A. From eq. (|A7J) it takes 
the form 

£v(t) = regular function — sinh (Or) sin (^t) • hi (u; r) , |argr| < n (27) 

for arbitrary complex r away from the cut. Evaluating eq. © for v < by closing the 
integration contour in the lower half t-plane one therefore encounters no singularities and 
the structure function vanishes identically. The logarithmic cut of the Vineyard function 
also shows up in the \ow-t expansion (see eqs. (jA18|) . (|A19p ) 



Zvttt) = MO) - -7- + -r— t 2 \n(iu Q t) 
Am Amn 



i 2 



,2 " 2 \ f^-\ .o/3 



- — J arctan I — I + 7O ^- - 7^ 



+ 0(t 3 ) (28) 



which has important consequences: first, we see that also the energy- weighted sum rule 
(EWSR) 

poo Qr 2 

j Q duuS(q,u) = z$;(0) = |- (29) 

is conserved whereas higher energy moments of the structure function diverge. This high- 
energy tail reflects, of course, the insufficient suppression of high frequencies in the simple 
model of Ohmic dissipation which may need modification for phenomenological applications 
in nuclear and quark physics. Second, it is even possible to determine the tail for large 
energy transfer analytically: as shown in Appendix B one obtains 

3 ^ M ~tJ_£ + ^ + fc). (30) 

Compared with eq. (|T§j) one sees that the suppression of the low- lying states by the square 
of the elastic form factor has disappeared and the asymptotic form (jSUj) does not depend 
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anymore on ujq or the oscillator parameter b - a property which roughly resembles the 
conjectured "quark-hadron duality" Q|. 

Another consequence of eq. (J28|) is that logaritmic corrections to ^/-scaling will persist 
even for large momentum transfer. This is because interaction times t ~ 1/q are probed in 



that limit 



28| and therefore q 2 t 2 In (iuj Q t) remains unbounded for q — > oo. 



IV. NUMERICAL RESULTS 

Sticking to pure Ohmic dissipation we next try to evaluate the structure function quan- 
titatively. Only for q — > (i.e. photoabsorption) the expansion (fT^j) into a sum over excited 
states is useful. For arbitrary momentum transfer the numerical problem is much harder as 
one has to calculate the inelastic structure function 

^etticfe v) = ^ exp[-2g 2 ^(0)] j ^ dt e wt { exp [2q 2 ^t) ] - 1 } (31) 

after subtraction of the elastic line (a 5-function) as a Fourier transform over an infinite 
interval. One may alleviate the numerical problem slighthly by expressing the inelastic 
structure function as a sine transform over the imaginary part of the characteristic function 
as demonstrated in ref. [12J 

S&Sedq,") = e"^ V - f°°dt sin(z4) (-Im) {exp [2q 2 fr(it)] } . (32) 

n Jo L L J > 

This holds since ^y(^t) vanishes as 1/t 2 in a sector of the complex t-plane which includes the 
lower half-plane (see Appendix A). Thus one may write a Cauchy integral representation for 
the (inelastic) characteristic function and express its real part in terms of the imaginary part. 
We have used the adaptive integration routine D01ASF from the NAG library together with 
the explicit representations (jA8|) . (jA13|) of the Vineyard function to perform the numerical 
evaluation of eq. (j22J)- 

Figs. Q and El show the results of the calculation for several momentum transfers and 
damping parameters. It is seen that the excitation of individual levels gradually moves into 
the broad structure of the quasi-elastic peak as the momentum transfer increases. The value 
7/cjo = 0.2 corresponds roughly to the one used in ref. 2^| where a parametrization of 
the response function was fitted to photoabsorption and electron scattering data in 12 C. 
A peak position of 22.7 MeV and a FWHM of (2a/1ii2 ■ 2.6 = 4.3) MeV was obtained for 
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FIG. 1: The structure function of a damped harmonic oscillator as a function of the energy transfer 
for momentum transfers qb^ = 0.5 (top) and qbo = 1 (bottom). Note the different scales in both 
plots. The solid curves are for a value of the Ohmic damping parameter 7/cJo = 0.2, the dashed one 
for 7/^0 = 0.3 and the dotted one for 7/cJo = 0.4. The undamped oscillator length and frequency 
are denoted by 60 and wo, respectively. 

the giant resonance. Although accounting for the coupling of this excitation to many-body 
states and the continuum this value of 7/u; gives too much structure from individual levels 
at higher q compared with typical experimental cross section. This may indicate that a 
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modification of the assumed Ohmic damping is needed for a description of medium-energy 
inclusive scattering data from nuclei. 





FIG. 2: Same as in Fig. ^but for qbo = 2 (top) and qbo = 4 (bottom). 

That the sum rules are well preserved can be seen in Tables |H] where the relative size of 
different contributions is presented for / y/u>o = 0.2 (other values of the damping parameter 
yield similar results). The numerical part was obtained by integrating the structure functions 
shown in Figs. ^ El with weight v n , n = 0, 1 from v = to v — v max in steps of Au by 
means of a Simpson rule. The NEWSR, of course, also gets a contribution exp(— q 2 b 2 /2) 
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from the elastic line. Also listed are the asymptotic contributions obtained from eq. (|3U|) 
by integrating from z/ max to oo. It is seen that the asymptotic contributions substantially 
improve the convergence to the correct values except for a few cases at high momentum 
where v max is not large enough. 

Due to its weighting the EWSR is, of course, not as well fulfilled as the NEWSR, but the 
agreement is very satisfactory and sufficient to demonstrate the consistency of all parts of 
the calculation: elastic line, inelastic excitations and asymptotic behaviour of the structure 
function. 
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num. 
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0.01 
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0.00031 
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0.62471 
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0.99193 




5. 
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0.00127 


0.00025 


0.99957 




7. 






0.37448 


0.00065 


0.00009 


0.99993 


2.0 


5. 


0.02 
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0.80314 


0.00509 
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0.96460 
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0.84239 


0.00200 


0.00099 


0.99768 
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0.84532 


0.00127 


0.00051 


0.99940 


4.0 


15. 


0.05 


0.00054 


0.97227 


0.00226 


0.00241 


0.97749 




20. 






0.99455 


0.00127 


0.00102 


0.99738 




25. 






0.99744 


0.00082 


0.00052 


0.99932 



TABLE I: Test of the non-energy-weighted (NEWSR) sum rule © for different momemtum trans- 
fers q and upper limits v ma x of the numerical integration performed in steps of Aza The parameter 
for Ohmic damping is taken as 7/wo = 0.2. The columns labelled "elast." and "num." denote 
the contribution of the elastic line and the result from numerical integration up to v = v maX) re- 
spectively. The columns "asy." list the leading-order (LO) and next-to-leading (NLO) asymptotic 
contributions from v max to infinity. In places with no entry the previous value applies. 

Although the NAG routine D01ASF does an excellent job in evaluating the oscillating 
integral over the characteristic function one may ask whether it is possible to introduce ad- 
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ditional damping by deforming the integration contour in eq. (J3~T|) in the upper-half t-plane. 
Indeed, running along both sides of the cut would eliminate the oscillating exponential factor 
altogether. However, this is not possible since as shown in Appendix A there are Stokes lines 
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0.9962 
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0.0085 


0.0068 
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20. 




0.9819 


0.0064 


0.0038 


0.9921 
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0.9898 


0.0051 
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TABLE II: Same as in Table [I] but for the energy-weighted (EWSR) sum rule (|2U)) divided by 
q 2 /{2m). Note that no elastic contribution exists in this case. The Ohmic parameter is again 
7/wo = 0.2. 



in the upper-half t-plane which separate the power-like decrease of the Vineyard function 
from an exponential increase which would overwhelm the exponential damping from 
the factor exp(ii/t). The optimal damping which is achievable without contributions from 
the arcs at infinity is a rotation by an angle 



<Pi 



arctan 



\2Q 



(33) 



which is depicted in Fig. 01 We have evaluated 

d.h.o. („ .A _ e ~^q 2 b 2 1 



5, 



inelastic 
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( dx exp ivxe w { exp 2q 2 £ v (r = ixe w \ — 1 1 

«/ 



(34) 
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FIG. 3: Integration contours in the complex t-plane for the numerical evaluation of the structure 
from the characteristic function: original path along the real axis (solid line) and path rotated by 
an angle ip into the upper half-plane (double line). The arcs at infinity (dashed curves) do not give 
a contribution as long as ip < ip\ where ip\ is defined in eq. ()33|) . The logarithmic cut along the 
positive imaginary axis is shown as hatched strip. 

by standard Gaussian integration after mapping the infinite interval to a finite one. As 
expected the damping factor exp(— vxsvn tp) is most beneficial for large energy transfer 
whereas a greater number of integration points is needed to avoid a negative numerical result 
for the structure function at small v. Thus the contour rotation method is complementary 
to an explicit summation over excited lines as given in eq. ([190. 

V. SUMMARY 

We have shown that the Caldeira-Leggett model of the damped harmonic oscillator also 
provides a consistent model for inclusive processes where it accounts for the coupling of a 
single particle to more complicated states and additional degrees of freedom. With only 
one additional (damping) parameter in the simplest version with purely Ohmic dissipation 
a highly non-trivial structure function has been obtained which conserves the sum rules. 
Its characteristic function has been given in closed analytic form with a clear separation 
between momentum transfer and "time", the variable conjugate to the energy transfer in 
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the process. This allowed us to concentrate on the study of a function of a single variable 
- the Vineyard function - and its analytic properties which determine the dynamics. As 
an extension to the many-body case is straight-forward further applications in nuclear or 
hadronic physics seem to be possible. 



Acknowledgments 

This work originated from a course "Path Integrals in Quantum Physics" given in the 
winter semester 2002/3 at ETH Zurich. I am grateful to all students who forced me to 
think about problems which I hadn't considered before. Many thanks also to Gerhard Baur 
(Julich) who pointed out relevant work on damping in nuclear reactions and multi-phonon 
giant resonances. 



APPENDIX A: THE VINEYARD FUNCTION FOR THE HARMONIC OSCIL- 
LATOR WITH OHMIC DAMPING 

Here we derive and collect a few properties of the Vineyard function defined by eqs. (fTol 
IT5|) in euclidean time r > 

1 roc cos(Er) 



£y( T ) = / dE— - 1 — ~ . (Al) 

^ VK ' 2mvrio E 2 + ^E + Uq K ' 

One may decompose the cosine function into exponentials and deform the integration path 
such that it runs along the imaginary axes. There is no contribution from poles of the 
integrand (which are all in the left-hand .E-plane) and one therefore obtains (first for r > 0) 

ry POO TTJ 

= S^i iE W^WT^ e ' (Rer - 0) - (A2) 

This is suitable for analytic continuation t —> it since the integral also converges for Rer > 0. 
Using partial fractions one gets 

1 ^ roo i 



£ V ( T ) = y rs / dE- — e 
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-Et 



]T rse E ^ T E 1 {E r<a r) , r>0 



r,s=±l 



with 



Here 



rft + is 



7 



EAz) 



n=l 



-z 



nn\ 



r, s — ±1 



=: —7s — In z + Ein(^) 



(A3) 

(A4) 
(A5) 



is the standard exponential integral and je — 0.57721566.. Euler's number. A careful 
evaluation of the arguments of the logarithm for r > then gives 
1 



£v(r) 



8rrmQ 



E 

r,s=±l 



f /2ft\ 


71 


< arctan 1 — 1 


— r — \- irs 




2 



7b + ln(a;or) - Ein (E rjS r) 



(A6) 

which defines the Vineyard function for arbitrary complex r with |argr| < it. Since Ein(z) 
is an entire function (ref. j^], chapter 5.1, footnote 3) one sees that 



£v( r ) = regular function — 



— sinh (ftr) sin f—r ) ■ In (ujqt) 

2mvrft K 1 V2 1 



(A7) 



which allows to determine the discontinuity across the cut. 

For purely imaginary arguments r = it ,t real, eq. (|A6|) can be written as 



r(t) 



1 / 7 

— exp — iVtt \t\ 

4mft F V 2 1 



r(t) 



Im 



f^) El {z 1 {t))-{ 1 ^- 1 )] ,z 7 (t) 



a'ft + ~) 1/1 



(A8) 
(A9) 



AranQ, 

The first term corresponds to the narrow-width approximation (J25|) whereas the last term 
(which is real, even in t and vanishing for 7 = 0) corrects for its deficiencies. 
Another representation of the remainder function r(t) is obtained by using the identity 



E 



(E 2 — Uq) 2 + 7 2 _E 2 



1 

4ft 



1 



1 



(A10) 



(£-ft) 2 + 7 2 /4 (£ + ft) 2 + 7 2 /4 
which shows that the weight function p(E) can also be considered as an (anti-)symmetrized 
Breit-Wigner distribution around E = ft . Inserting this into eq. (|A2|) one obtains for 
t — it, t real 

7 



8m7rft 



7 

8m7rft 



+°° exp(-iEt) 
-00 (£-ft) 2 + 7 2 /4 



2tt 

7 



exp 



-int 



7, 



expj-zgg 

-00 (E - ft) 2 + 7 2 /4 

°° exp(-zEt) 1 

(£ + ft) 2 + 7 2 /4 J 

f°° cos(Et) 1 

io (£ + ft) 2 + 7 2 /4 J 



(All) 
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Hence 

""'^r ^i^ ' 1 - 1 ' < ai2) 

Again one may distort the integration path such that it runs along the imaginary axis. 
Decomposing the cosine function into exponentials and realizing that the integrand in eq. 
(|A12|) has only poles at E — —Q ±27/2, one then obtains a representation 

r(t) = - 1 — / dE— =- — — e~ m , treal (A13) 

v ; 2rmr Jo (E 2 - ul) 2 + 4tt 2 E 2 ' v ' 

which is very well suited for numerical evaluation. Indeed, we have checked the routine which 
calculates the Vineyard function £y(it) based on on the exponential integral representation 
2U) by a direct Gaussian integration of eq. (jA13|) and found a relative deviation 



< 2 • 1(T 6 (A14) 



for all real t and 'j/ujq < 0.4. In this comparison the complex exponential integral Ei(z) was 
calculated by the rational approximations with n — 10 terms given in ref. |3J|: Table 64.4 
was used f or \z \ < 9 and Table 64.5 for \z\ > 9 and checked against values listed in Table 
5.6 of ref. |3i|. Fig. EJshows £v( r ) (relative to the undamped case blexp(— lu t)/A ) and 
Fig. El the remainder function r(t) for the chosen values of the damping parameter 7. 
Note that one has a simple result for the imaginary part (which is odd in t) 

lm£ v (it) = -— — sin (fit) e" 7l ' l/2 (A15) 

but that the real part is more involved. This can also be seen if we use the standard 
representation of the step function 

e(E) = -— f + °° ds *M-isE) (A16) 

v 1 2ni J-co s + ie K J 

to extend the integration range in eq. (jA2|) to — 00 . Using the identity (|A10|I simple 
manipulations then give 

tvfr) = j 1 ^ rds^^-e-^ , treal (A17) 
AmTT\l J -00 s — t + te 

which shows that the real part is determined by a principal-value integral whereas the 
imaginary part is given by eq. (|A15JI . 
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FIG. 4: The Vineyard function £v(t) normalized to the undamped case as function of the euclidean 
time r. The solid curve is for a value of the Ohmic damping parameter 7/wo = 0.2, the dashed 
one for 7/uo = 0.3 and the dotted one for 7/^0 = 0.4. 

Finally we consider the behaviour of the Vineyard function for small and large r. From 
eq. (|A6|) one immediately obtains 



£v(t) = E a fc rfe + E b ^r 2k ln(w r) 

fc=0 fc=l 



(A18) 



with 



a 



d-2 



1 (2& 

— arctan — 

2mn{l \ 7 



a 1 



Am 



1 



,-w3 \ f„ 2 7 2 \ 
jil [ - - 7^J + I il - — I arctan I — — 



a 3 



<^ _ ^\A19) 
48m V if 1 



7 



7 



4m7r 



24m7r 
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FIG. 5: The correction term r(t) to the narrow-width approximation (|A8|) for the Vineyard function 
as function of the time t. Note the logaritmic scale for r(t). 

Another possibility is to use the differential eqaution 

&(t) - (2ul - 7 2 ) #(r) + uj^ v (t) = JL- 1 (A20) 
v > Amu T z 



with the appropriate boundary conditions which follows directly from eq. (jA2l) . 

The asymptotic behaviour for arbitrary r is more involved. Eq. (jA2|) may be used in the 
right-hand r-plane: one simply has to expand p(E) for small E and integrate term by term 
to find 

Zv(t) ™ 7^^ + —* ^ 2 -t) V--- Re ->°- (A21) 
ZrrnTujQ t z mirujQ \ 4 I r 4 

This is consistent with the result from the differential equation (jA20|) for large r if one 
considers the derivatives as corrections. To find the asymptotic behaviour in the left-hand 
r-plane we may use the explicit representation (|A6|) and re-introduce the exponential integral 
as this function has a simple asymptotic behaviour (see, e.g. eq. 5.1.51 in ref. jscl ]) 



Ei(z) 



1! 2! 

- + - 

z z z 



\aigz\ < 3tt/2 . (A22) 
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However, care is needed when replacing Ein(z) by the exponential integral since the correct 
addition theorem of the logarithm with complex arguments 

In (aft) = ln(a) + ln(6) + 2ni [e(-Ima)6(-Im 6)6(Im (ab)) - 6(Ima)e(Im6)e(-Im (ab)) 

(A23) 

has to be used for \n(E r ^ s T). Consequently 

£ v (t) = £ re E ^\27fQ(slmr)Q(-slm(E r>s T))+tsE 1 (E r>s r) 



arerr < 7r 



(A24) 

Although eq. (|A24|) is less suited to display the analytic structure of the Vineyard function 
it allows to find the Stokes lines for the asymptotic behaviour. These are the rays in the 
complex r-plane which divide the power-like decrease of eq. (|A21|) from an exponential 
increase. Indeed, due to eq. (jA22|) the last term in the square brackets of eq. (jA24|) gives 
rise to the power-like decrease of £v( T ) but this will be overwhelmed by the exponential 
increase of the first term if 

e(Re(£ r , s T)) e(sImr)e(-sIm(£ r , s T)) ^ (A25) 

for any r, s = ±1. Writing 

E h i = VL + i^ = ujoe^ 1 , ipx = arctanQ^) , < tp x < | (A26) 

a straightforward analysis of the condition ()A25jl shows that an exponential increase only 
occurs for 

|argr| > | + ^ , (A27) 

i.e. inside a sector around the cut with opening angle ir/2 — (p 1 = arctan(2fi/7). This means 
that eq. (jA21|) holds for the wider range |arg r| < ir/2 + <pi . 



APPENDIX B: ASYMPTOTIC BEHAVIOUR OF THE STRUCTURE FUNC- 
TION 



Here we derive the asymptotic behaviour of the structure function when the energy trans- 
fer v becomes very large. This is done by standard asymptotic analysis: for example, one 
may apply eq. (30) in ref. to our eq. (|H2^). Then one obtains 



S(q,u) 



1 



1 



- Im $ ? (0) + — Im <(0) - -= Im $J K (0) + 



7 q | 



(Bl) 
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since lm$ g (0) = 0. Here we have used the representation of the characteristic function 
in terms of the Vineyard function £y(it) and the low-t expansion of the latter. Note that 
the leading contribution comes from the logarithmic term in eq. (|28|) which produces an 
imaginary part for £y(0). 

However, higher-order terms cannot be calculated by means of eq. (jBlJ) since lm$^ y (0) 
does not exist. This shows that the next-to- leading asymptotic term is not falling off like 
1/zA To determine this term we use the exponential representation (jUJ), well-known results 
for the Fourier transform of generalized functions and the low-t behaviour of the Vineyard 
function. Writing 



2g 2 [£v(i*)-fr(0)l = : f(t) + g(t)\n(zuj t) = f(t)+g(t) 
with 



.7T 



In (u \t\) -M-sgnt 



/(*) 



AID. V 



2m 



t^o 7 2 

g(t) — ► * 

ii 2m 



O It 



(B2) 



(B3) 



one simply gets by expanding the exponential 



l + f(t)+g(t) 
1 



7T 



In (w |*|) + i-sgnt 



+^f 2 (t) + f(t)g(t) 



71 



In (£^ |t|) + i-sgnt 



(B4) 



+\g\t) 



7T 



In (uj \t\) + i-sgnt 



+ 



As the Fourier transform of powers of t gives derivatives of ^-functions we see that all regular 
terms do not contribute to the asymptotic behaviour for large v. The contribution of the 
non-analytic terms can be taken from Table 1 of ref. 32\ (setting y = —vj(2ix) ) 



+oo 



iut 



oo 
)-oo 



dt t n sgn t e 
dtt n In It I e ivt 



in+l 



m 

m- — , sgnz/ . 



(B5) 
(B6) 



One then realizes that the leading contribution in the asymptotic expansion arises from the 
last term in the first line of eq. (jB4|) 



r,/ \ i/— >oo 1 7<? 



2ix 2mn 



dtt 2 



71 



In (u \t\) -M-sgnt 



7 q 



— h > ( B7 ) 



mix u° 



in agreement with eq. IjBlj) . The subleading contribution stems from the last term in the 
second line 
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whereas one can show that the last line gives a contribution of order (lnz/)/iA The asymp- 
totic expansion makes sense if the subleading term is much smaller than the leading one 
which requires v ^> q 2 /(2m), i.e. excitation energies much larger than the maximum of the 
quasielastic peak. 
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